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Bloch oscillations of a soliton in a molecular chain 
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Abstract. The paper presents results of numerical experiments simulating Bloch oscillations of solitons in 
a deformable molecular chain in a constant electric field. By the example of a homogeneous polynucleotide 
chain it is shown that the system under consideration can demonstrate complicated dynamical regimes 
when at the field intensities less than a certain critical value, a soliton as a whole exhibits oscillations, 
while at the field intensities exceeding the threshold, a soliton turns to a breather which oscillates. It is 
shown that the motion of a charge in a deformable chain is infinite as contrasted to that in a rigid chain. 
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It is well known that an electron occurring in an ideal 
rigid periodic molecular chain or in a solid state superlat- 
tice exhibits Bloch oscillations in response to a constant 
electric field [l]-[5]. In an external time-periodic field, mo- 
tion of a charge along a rigid chain can be both infinite and 
finite (dynamical localization) [B]-[TU]. In a deformable 
crystal chain the role of an external field is played by os- 
cillations of the lattice nodes which can be presented as 
superposition of plane travelling waves, or phonons. In this 
case the motion of an electron along the chain is thought 
to be infinite since the electron scatters on phonons and 
Bloch oscillations do not take place [TT| . 

It is common knowledge that in quasi-one-dimensional 
molecular chains interaction of an electron with lattice os- 
cillations is not weak. Therefore we cannot safely assume 
that the electron wave function goes off phase (in view of 
scattering of the electron on phonons) and Bloch oscilla- 
tions fail. 

To clear up this point we consider the case when a 
charge placed in a molecular chain transits to a soliton 
state as a result of interaction with lattice oscillations. 
This occurs, for example, in homogeneous polynucleotide 
chains where the charge motion is described by Holstein 
Hamiltonian in which each site presents a nucleotide pair 
considered as a harmonical oscillator [T21-I14I: 
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Hh = i^ ^(a+a„_i + a+a„+i) +^ anO^an , (1) 

n— 1 n— 1 
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Here Hh - is a Hamiltonian of a charged particle, a^,an 
- are operators of creation and annihilation of the charge 
on site n, i^ - is the matrix element of the transition from 
the n - th site to the n ± 1 - site, a„ - is the energy of 
the particle at the n - th site, Tiujb = eSa, where E - is 
the intensity of the electric field, e - is the electron charge, 
a - is the distance between neighboring bases. Tk - is an 
operator of the kinetic energy of sites. Up - is the potential 
energy of sites, P„ - is an impulse operator canonically 
conjugated to the displacement g„, M - is the effective 
mass of the site, fc - is an elastic constant, a' - is the 
particle-site displacement coupling constant. 

We can pass on to semi classical description of the 
wave function of the system \^{t)) as an expansion over 
coherent states: 

N 

Wit)) =^&„(i)a+exp{~ J^[/3,(i)P,- - 7r,(i)g,] }|0), (2) 

where |0) - is the vacuum wave function and the quantities 
Pj{t) and 7rj(i) satisfy the relations: 

{^{t)\qr,W{t)) ^ f3,,{t), {9{t)\P^\^it)) = 7r„(i). (3) 

Dynamical equations for the quantities 6„(i) and /3n(i) 
resulting from (HJ - ([21) have the form : 
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Ctnbn + v{bn-l +6„+l), 



M/3„ = -7/3„-fc/3„-a'|6„|^ 



(4) 
(5) 
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Equations ^ are Schrodinger equations where 6„ is 
the amplitude of the particle localization at the n - th site. 
Equations ^ are classical motion equations describing 
dynamics of nucleotide pairs with regard for dissipation. 
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Fig. 1. Oscillatory motions of a soliton for some values of 
parameter k(k = 1, 2, 3, 4) at the electric field intensity _E = 0.1. 
The length of the homogeneous nucleotide chain is A'' = 201, 
t ^ t/r, r = 10-^*sec, S' = 0.006, Zj = 0.01, 77 = 1.276. 

(2' =ij'T, Lj' =7/M = 6- 10"sec-^ 'JT, I^ = Lor). 



where 7 is friction coefficient. We believe that a semiclas- 
sical description in which motion of a charge along a chain 
is described by quantum motion equations Q and motion 
of individual nucleotides is presented by classical motion 
equations ([5]) is valid in view of a large nucleotide mass 
(« 300 proton mass). 

In the case of a rigid chain, when a' = the solution 
of the system (g]), © wih be [11], [H]: 
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Jn{x) - is Bessel function of the first kind. Solution ^ 
corresponds to Bloch oscillations of a particle in the chain 
affected by an electric field for which the particle's centre 
mass: 
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X{t) = Y. K{t)\^na, 
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demonstrates periodic oscillations at the frequency of ujb '■ 
2av 



X{t) = X{0) + - — \So\ (cos 00 - cosiiOBt + 0o)), 
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where a is the distance between neighboring nucleotides, 
which for DNA is equal to SAA. 

For a' ^ in the absence of an electric field, a station- 
ary solution of equations ^ , dS]) corresponds to a localized 
state of a soliton type. To study the evolution of a soli- 
ton state in an electric field we will use an initial charge 
density distribution such that: 



|&„(0)| = — cosh ( ^^ 
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Initial values of x^ and y" (6„ — Xn + iyn) for i^ > have 
the form: 

4 = |6„(0)|(_l)yy2, y,", = |6„(0)|(-ir+yy2, (10) 



which corresponds to the gr^ 
absence of an electric field 

Fig.l shows the results 
Q, dS]) for some values of 
the intensity of the charge 
lations at the electric field 
^ = WT = 0.01, ?7 = 1.276. 
w and 77 are the same as in 
In dimensional units these 



ound state of a particle in the 

mm- 

of the solution of equations 
parameter k, responsible for 
interaction with lattice oscil- 
intensity E = Eear/Ti = 0.1, 
Here the values of parameters 
work [T^, and r = lO^^^'sec. 
parameter values correspond 



V.D.Lakhno, A.N.Korshunova: Bloch oscillations of a soliton in a molecular chain 



x<xy 



JCCf) 




1 15D 



1000 



Fig. 2. Transition of a soliton into a breather for k = 4. 
(Fig.l d) at large times.) 



io £ = 1.94 • Kf'V/cm, u = yJk/M = W^sec-^, 
V — 0.084el/. The parameter of electron-phonon strength 
K = 4, which in dimensional units corresponds to a' = 
Q.l?>eV/A is the same as in [12]. This value is close to 
that used by other authors (in [TB] a' was found to be 
a' w0.23el//i). 

It is seen from FiglT] a) that in the presence of an elec- 
tric field, a soliton executes periodical motions, coming 
back to the point where the soliton center mass initially 
occurred. This oscillatory motion corresponds to Bloch 
oscillations with a period of T = 2tt/ujb- The total am- 
plitude of the oscillations L is close to that determined 
from the solution of linear problem ^ and is written as: 
AWa/E, where AW = A-q stands for the width of the 
conductivity band equal to AWr/h in dimension form. 
For the parameter values presented above L « 51a (with 
the soliton characteristic size « 10a). 

FiglT] b)-d) shows evolution of the dynamical behavior 
of a soliton at the initial stages of the motion as parameter 
K increases. After a lapse of time Bloch oscillations restore 
(restored Bloch oscillations are not given in Fig|l]b)-d)). 

In the case of strong electric fields presented in FigH] 
a soliton executing Bloch oscillations with time turns to a 
breather oscillating at Bloch frequency (Fig|2]). At rather 
large values of n, a breather can arise from the initial soli- 
ton state immediately, i.e. by-passing the phase of Bloch 
oscillations as a whole. 

Without going into details of nonstationary regimes 
of the particle motion in the cases under consideration we 
will restrict ourselves to mere qualitative description of the 
picture. It has been observed that the case of a deformable 
chain (a' ^ 0) differs qualitatively from the limiting case 
of a rigid chain (a' = 0) in that at finite a' the quantity 
X{t) given by ([7]) grows infinitely at t — t' oo (Fig|3]). This 
result could have been guessed from the already mentioned 
analogy between the influence of a periodic external elec- 
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Fig. 3. Function X{t) for various values of uj'. 



trie field on a particle and oscillations of phonons. Quite 
nontrivial is the finding that under this infiuence, in the 
case of a strong particle-phonons interaction, i.e. when a 
soliton is formed, Bloch oscillations of the particle persist 
in the electric field as oscillations of a soliton as a whole 
or a breather, depending on the system parameters. 

In conclusion it may be said that this picture of the 
charge motion in a deformable molecular chain in a con- 
stant electric field at zero temperature T = seems to be 
rather general: a positive charge introduced in the chain 
will move along the field executing Bloch oscillations. At 
finite temperatures a soliton or breather state will break 
thus leading to failure of Bloch oscillations. In this case 
motion of the charge over the chain will be infinite along 
the lines of the field and have an ordinary band character. 
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